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ABSTRACT 

We present near-simultaneous Chandra/HST observations of the very faint ( L x < 10 36 erg s -1 ) X-ray tran¬ 
sient source M15 X-3, as well as unpublished archival Chandra observations of M15 X-3. The Chandra 
observations constrain the luminosity of M15 X-3 to be < 10 34 erg s -1 in all observed epochs. The X-ray spec¬ 
trum shows evidence of curvature, and prefers a fit to a broken power-law with break energy Ebreak = 2.7 

keV, and power law indices of Ti = 1.3^2 an d ^ = 1.9^ 2 over a single power law. We fit our new F438W 
(B), F606W (broad V), and F814W (/) HST data on the blue optical counterpart with a model for an accretion 
disk and a metal-poor main sequence star. From this fit, we determine the companion to be consistent with a 
main sequence star of mass 0.44friyjjg 3 , M 0 in a ^4-hour orbit. X-ray irradiation of the companion is likely to 
be a factor in the optical emission from the system, which permits the companion to be smaller than calculated 
above, but larger than 0.15 Mq at the 3 ct confidence level. M15 X-3 seems to be inconsistent with all suggested 
hypotheses explaining very faint transient behavior, except for magnetospherically inhibited accretion. 

Subject headings: accretion, globular clusters: individual (M15), stars: neutron, techniques: photometric, X- 
rays: binaries. X-rays: individual (M15 X-3). 


1. INTRODUCTION 

The Galactic globular clusters (GCs) are host to many ex¬ 
otic o bjects, due to the high central stella r density at their 
cores dClarkl 1 1 9751: iVerbunt & Lewinl120061) . One such class 
of objects are the low-mass X-ray binaries (LMXBs), sys¬ 
tems that consist of a compact object accreting material from 
a companion overflowing its Roche lobe. Many LMXBs ex¬ 
hibit transient behavior, where their luminosity can change by 
a factor of 100 or more. This behavior is generally believed 
to be_tied to_ thermal-viscous instability in the accretion disk 
(lLasotall2001l) . 

One broad category that is used to classify transient X- 
ray sources is their peak X-ray luminosity. Sources that 
have a peak luminosity exceeding 10 37 erg s~ are clas¬ 
sified as bright sources. With the deployment of sensi¬ 
tive telescopes including Chandra, XMM-Newton, and Swift 
in the past fifteen years, the class of very faint LMXBs, 
whose peak X-ray luminosities lie be l ow 10 36 erg s -1 , have 
been discovered dMuno et ~akl l2005bl IWiin ands et alJ 12006b 

lHands et all 12004 IDegenaar & Wiinandsl 120091) . These ob¬ 

jects present a puzzle for binary evolution, as their luminos¬ 
ity implies a time-averaged accretion rate that can reach as 
low as 10~ 13 Mgyr” 1 or less ([King & Wiinandsll2006l) . Sev¬ 
eral explanations have been offered for these low accretion 
rates, including accretion from a planetary or brown dwarf 
comp anion, or by a prim ordial intermediate-mass black hole 
( King & Wiinandsl2006l) . accretion in an ultracompact system 
(lin’t Zand et al.l 120051). inefficient accre t ion due to magne- 
tospheric inhibition dHeinke et all 120091 : ID’ Angelo & Spruitl 
120121 [Heinke et alj|20l4. or accretion from a companion in 
the period gap dMaccarone & Patrunoll20l3h . Some VFXTs 
can be explained as an intrinsically bright object viewed at an 

1 Department of Physics, University of Alberta, CCIS 4-181, Edmonton, 
AB T6G 2E1, Canada; ramason@ualberta.ca 

2 Max-Planck-Institute fur Radioastronomie, Auf dem Hiigel 69, 53121, 
Bonn, Germany 

3 Astronomy Dept.. Indiana University, 727 East 3rd St., Bloomington, 
IN 47405 


unfavour able inclination, but this is not possible f or all such 
systems dMuno et alJl2005al IWiinands et al.ll2006l) . Observa¬ 
tions of VFXTs in optical and infrared can constrain the na¬ 
ture of VFXT systems by identifying the companion in the 
system, so that the unexplained accretion regimes of these 
sources can be better understood. 

Recent X-ray monitoring campaigns of the Galactic Cen¬ 
ter have measured the long-term behavior of a num¬ 
ber of VFXTs, identifying roughly as many objects that 
show only VFXT out bursts as tho se that show "nor¬ 
mal", bright o utbursts (IDegenaar & Wiinandsl 120091 120101 
IDegenaar et al.l 120121) . Short, low-luminosity VFXT out¬ 
bursts are also seen from transients showing bright outbursts 
(e.g.. IDegenaar & Wiinandil2009l) . A few systems are seen 
to maintain quasi-persistent (for at least several continuous 
years) X-ray e mission, well below 10 36 erg s -1 , incl uding 
NGC 6652 B dHeinke et al l |200H IStacev et al.l [20l2h, A X 
J1754 . 2-2754 (fSa kano et ak J20021 IChelovekov & Grebenevl 
[20071 IDegenaar et all I201HL 1RXS J171824.2-402934 


(in’t Zand et al 


dDel Santo et al.l 


2005, 


_ X MMU J174716.1-2 81048 

_ 1200% . M 15 X-3 dHeinke et all[20091) . and 

1RXH J 173524.4-354013 (IDegenaar et al.ll2010l) . These sys¬ 
tems are too faint to keep normal-sized (i.e., few ho urs o r¬ 
bital period) accretion disks irradiated, leading lln’t Zand et all 
(120051) to propose that the orbital separation must be ex¬ 
tremely compact (orbital period of 7 minutes in one case); 


however, such a short orbital period would produce a 
very high mass transfer rate (e.g.- lPelove & Bildstenlf2003 j). 
leading to a contradiction. Additionally, the measure- 
ment_of _strong Ha emission lines in at least one source 
(IDegenaar et al.ll2010l) shows a hydrogen-rich donor star, rul¬ 
ing out an ultracompact system in that case. 

The X-ray spectra of Galactic Center VFXTs have typically 
been well-described by a p ower-law of photon index rang¬ 
ing from ~1 - 3 ( iDegenaar & Wiinandsl2010l: IDegenaar et all 
120121) . However, the spectra have typically not been of 
high quality due to the high absorption, the small effec¬ 
tive area of Swift/XRT, or the short exposures in Chandra 
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or XMM-Newton monitoring observations. Three of the 
quasi-persistent VFXTs listed above have long XMM-Newton 
exposures at luminosities of 10 34 to 10 35 erg s~ 1 , show¬ 
ing that all have rather soft spectra, with photon indices 
T=2.3-2.5, but that the objects fainter than 10 35 erg s -1 re¬ 
quire a blackbody-like comp onent in addition to a power-law 
(1 Armas Padilla et all 12013 al). Another transient VFXT, IGR 
J17494-3030, observed with XMM-Newton at L x ~ 8 x 10 34 
erg s -1 , also requires a blackbody plus power-law in its spec¬ 
tral fit (lArmas Padilla et ahll2013bl) . 

The high extinction towards most VFXTs has hampered 
deep searches for optical or infrared counterparts. However, 
a few likely or certain optical/infrared counterparts have been 
identified. iDeutsch et al.l ( 1998 1 identified a relatively blue 
obje ct as a likely LMXB counterpart in NGC 6652, which 
IHeinke et al.l (12001 1) showed to be the counterpart to NGC 
6652 B; this object is extr emely variable in t h e optical and X- 
ray on minute timescales (lEngel et alil2012 1. iDegenaar et al.l 
(1201 Oh found an optical counterpart for 1RXH J173524.4- 
354013, which shows a strong hydrogen emission line and 
magnit udes characteristic of a low-mass main sequence star. 
iKaur et aTI (1201 2b identified a likely variable near-infrared 
counterpart to XMMU J174716.1-281048, indicating a low- 
mass star. The likely VFXT, and black hole candidate. Swift 
J 1357.2-0933 h as a low extinction, and a potential M4 star 
companion dRau et al.l 1 2011 1 although this potential com¬ 
panion was not detected when the system was quiescent 
(IShahbaz et al.ll20l3l) . 

M15 X-3 was discovered from archival Chandra High Res¬ 
olution Camera (HRC-I) and High Energy Trans missio n Grat¬ 
ing Spectrometer (HETGS) observations, by IHeinke et al.l 
( 2009 1. This source, located roughly 21" from the center of 
the globular cluster M15, has been observed at an Lx of either 
ss 10 34 erg s” 1 (in its bright state) or « 2 x 10 31 erg s” 1 (in 
quiescence). In zeroth-order HETGS observations M15 X-3 
has demonstrated a relatively hard, power law-like spectrum 
(r = 1.5) in it s bri ght state. It is likely to be a softer source 
in quiescence (IHeinke et al.ll2009[) . though a spectrum of the 
faint state of M15 X-3 would require a long (> 50 ks) observa¬ 
tion with Chandra triggered when M15 X-3 is known to be in 
quiescence. 'Heinke et al.l ( 2009 ) identified the optical coun¬ 
terpart to M15 X-3 as a likely main-sequence star, though the 
photometry had rather low signal-to-noise ratios. Throughout 
this analysis, we assume that M15 is at a distance of 10.3 kpc, 
with a co lumn density of 4.6 x 10 20 cm -2 for consistency with 
IHeinke et al.l (120091) . an E (B—V) of 0.10 and a cluster metal- 
licity [Fe/H] of -2.37 dJanulisI fl99l Ivan den Bosch et ahl 

l2006h . 

2. X-RAY DATA ANALYSIS 

2.1. ACIS-S Data Reduction 

For this analysis, we reduced and analyzed the 2012 Ad¬ 
vanced CCD Imaging Spectrometer (ACIS-S) Chandra obser¬ 
vation of M15, taken in subarray mode with an effective frame 
time of 0.9 s. We also reduced and analyzed three unpublished 
archival observations from the Chandra data archive (see 
Table []}• These observations were reduced using CALDB 
version 4.5.9 and the August 2012 time-dependent gain file 
(acisD2012-08-01t_gainN0006.fits). All four observations 
were taken using ACIS-S in faint mode with no grating. Each 
observation was reprocessed with the chandra_repro script, 
to generate a new level 2 events file. Using CIAO 4.6, the 
spectrum of M15 X-3 was obtained for each observation via 


Table 1 

Summary of Chandra Observations of M15 X-3 


ObsID 

Date 

Exposure 

Instrument 

L x 



[ks] 


[10 33 erg s _1 ] 

9584 

2007-09-05 

2.15 

HRC-I 

6 -2 

11029 

2009-08-26 

34.18 

ACIS-S 

o t+0.4 
-0.4 

11886 

2009-08-28 

13.62 

ACIS-S 

Q n+0.7 
yy -0J 

11030 

2009-09-23 

49.22 

ACIS-S 

q 2+0.4 

y ' z -0.4 

13420 

2011-05-30 

1.45 

HRC-I 

5+ 2 

3 -\ 

13710 

2012-09-18 

4.88 

ACIS-S 

I oil 


the specextract script. In all observations, M15 X-3 can be 
clearly detected as a source 21" from the core of M15 (see 
Figure Q}. Two other bright X-ray s ources, the LMXBs AC- 
211 and M15 X-2 (White & Angelinill2001h . are also nearby 
in the Chandra frame. The source region used was a 2.4" cir¬ 
cular region, surrounded by an annulus with 6" radius that 
provided a representative background region. Background 
regions of varying sizes were tested, however the choice of 
background region did not affect our results. The spectra were 
grouped to give a minimum of 25 counts per bin in the 0.5- 
10 keV range, so that Gaussian statistics were valid. Each 
spectrum was then fit using XSPEC 12.8.0. Initial fitting in¬ 
dicated that M15 X-3 was in a consistent spectral state across 
all four observations, motivating the choice to fit the spectra 
simultaneously, first with model parameters tied and then un¬ 
tied from each other. To perform a check for variability in 
M15 X-3, background subtracted light curves were generated 
using the dmextract command. Testing for variability using 
the CIAO tool glva nfl which looks for devia tions among op¬ 
timally sized bins ([Gregory & Loreddl 1 9921) . indicated evi¬ 
dence of variability within all 4 observations (probability of 
nonvariability always less than 10 8 ), as p re vious ly noted in 
other observations of M15 X-3 (IHeinke et al.1120091) . 

2.2. HRC-I Data Reduction 

There are five HRC-I observations of M15 X-3 in the 
Chandra archive. The first three, obtained in 2001, found 
M15 X-3 in a quiescent state, barely detectable when using 
special processing (IHeinke et alJ 1200 9). and we do not re¬ 
analyze them here. We analyze the 2007 and 2011 HRC-I 
observations, which catch M15 X-3 in a relatively bright state 
(see details in TableQ}. We reprocessed the events files (with 
CALDB 4.5.9 and the August 2012 time-dependent gainfile) 
using the chandra_repro script, and extracted source counts 
from the position of M15 X-3, and background counts from 
an annular region. We infer the unabsorbed 0.5-10 keV lu¬ 
minosities using PIMMS and a power law with photon index 
T = 1.4 and N H = 4.6 x 10 2l) cm -2 . The short length of these 
observations (Table [[]) makes the data unsuitable for timing 
analyses. 


2.3. Ml5 X-3 Spectral Analysis 

Previously, M15 X-3 was found to be well-described using 
a power law with photon index 1.51 ± 0.14, absorbed through 
a column density fixed at the cluster value of 4.6 x 10 20 cm 2 

4 http://cxc.harvard.edu/ciao/ahelp/glvaiy.html 
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Figure 1. X-ray image of M15 (ObsID 11029). The inner circle (radius 2.4") 
indicates the source extraction region for M15 X-3, and the inner bound¬ 
ary of the background extraction region. The outer circle (radius 6") is the 
outer boundary of the background extraction region. Note the presence of the 
heavily piled up LMXB sources AC-211 (left) and M15 X-2, as well as the 
Chandra ACIS-S readout streak. M15 X-3 lies approximately 2 1" from the 
center of M15, well outside the < 1.5" core ISosin & Kind 19971) . 

dHeinke et al .1120091) . Very faint X-ray transients are often de¬ 
scribed with a hard powe r-law, or a hard power-law with a 
softer thermal component (lArmas Padilla et aT]l2013al) . Even 
for a relatively faint source such as M15 X-3, the count rate 
of « 0.06 counts s _I (or 0.18 per frame, for the three long 
archival ACIS-S observations), makes pileup a potential is¬ 
sue. Therefore, pileup is accounted for by convolving spectral 
models with the XSP EC model for pileup in Chandra obser¬ 
vations (lDavisll2d0fl) . The dependence of fits upon the grade 
correction parameter a was considered in the fitting process 
0 . 

For this analysis we fit all four observations simultaneously. 
These separate fits all suggested that M15 X-3 was in the same 
state with L A ps 8 x 10 33 erg s _1 (\ 2 = 211.73 for 188 degrees 
of freedom). This motivated the choice to fit all four ACIS 
observations of M15 X-3 simultaneously, with the power-law 
index of all four observations tied together. In this case, only 
the power law component normalization, which corresponded 
to the measured flux of the source, was allowed to be indepen¬ 
dent. Fitting to all of the four Chandra observations gives an 
acceptable fit (x 2 /dof = 214.13/191, implying a null hypoth¬ 
esis probability of 0.12, see Table[2]SPL 1) with an absorbed 
power law of photon index 1.42 ± 0.03. Allowing the indices 
to vary does not significantly improve the fit (an F-test gives 
an F-statistic of 0.71, and probability 0.55 for obtaining such 
an improvement by chance), suggesting that the power-law 
index does not change. The fit was not significantly affected 
by the value of the grade correction parameter a, which was 
unconstrained in the fit, due to the relatively small degree of 
pileup. We therefore report results with a fixed to 0.5, as em¬ 
pirically suggested ( IDavisI120011) , as none of the results pre¬ 
sented here qualitatively change for different choices of a. 

Examining the single power-law fit suggests the presence of 

5 N a gives the probability that the photon will be identified as a ’bad' grade 
and rejected, where N is the number of photons landing in the detection cell 
during one frame time lDavlsi200Th . 


curvature in the residuals above 1 keV. This curvature moti¬ 
vated the choice to attempt other fits. First, the column density 
Nh was allowed to vary, since high column density is a poten¬ 
tial source of curvature. All subsequent fits permit Nh to vary. 
In the case of the single power-law, this strongly improved the 
fit (x 2 /dof = 195.08/190; SPL 2 in Table 0 with a new col¬ 
umn density of 9j/ 2 x 10 20 cm -2 , though the residuals of the 
fit (Figure O still suggested curvature. Since other VFXTs 
are described acceptably by a s pectra composed of a power - 
law with a thermal component (lArmas Padilla et ahll2013al) , 
several thermal components (a low magn etic field n eutro n 
star hydrogen atmosphere, NSATMOS, iHeinke et al.l (i2009t): 
a mu lticolor disk blackbody, DISKBB, e.g. iMitsuda et all 
(119841) : and a single-temperature blackbody, BBODYRAD) 
were tested to see if they improved the fit. 

Adding an NSATMOS component (assuming emission 
from the entire neutron surface, i.e., normalization fixed to 
1.0; NSATMOS + PL in Table 0 does not substantially im¬ 
prove the fit (A\ 2 is only 0.04 for 1 fewer d.o.f), and the fit 
prefers a temperature of 0.051 qq 4 keV. The 90% confidence 
upper limit on the temperature permits this component to con¬ 
tribute only 2% (~2x 10 32 erg s _1 ) of the 0.5—10 keV unab¬ 
sorbed Lx, and the lower limit is pegged to the lower tempera¬ 
ture boundary of the model. Adding a disk blackbody compo¬ 
nent (DISKBB + PL in Table 0 does improve the fit (\ 2 /dof 
= 188.04/188, F-statistic of 3.52 and probability 3.2 x 10~ 2 
of this improvement being due to chance). However, the disk 
blackbody model can be physically excluded, since the in¬ 
ferred inner disk radius (R mncr cos 9 = 0.7 km) is much smaller 
than the neutron star radius, making an accretion disk origin 
untenable except at very high inclination. 

Addition of a blackbody model can also improve the fit 
(X 2 /dof = 187.52/188, F-statistic of 3.79 and probability 
2.4 x 10 2 of this improvement compared to the variable Nh 
single power-law model being due to chance). The best-fit 
blackbody parameters are kT=Q.7+QQ7 keV, with an emitting 

radius R=0.2+q ] km (see Table [2] BB + PL). This could 
possibly be interpreted as a hot spot on the neutron star 
surface, where infalling material is channeled by the mag¬ 
netic field to the poles. However, the inferred emitting re¬ 
gion is smaller than is expected for polar caps on a neu¬ 
tron star surface spinning at millisecond periods (1.4-2.6 
km for 3-1 0 ms periods, using R pc = (2 ttR N s/(cP)Y^ 2 Rns', 
[Lvne & Graham-Smithl2006h . A longer spin period (e.g., 0.5 
seconds for 0.2 km) would predict a polar cap of the appro¬ 
priate size. 

Altern ative ly , it is well-known (e.g. IZavlin et al] 119961 ; 
iRai agonal & Romanil Il996i) that a hydrogen atmosphere 
model gives inferred radii a few times larger than black¬ 
body fits to the same data. Indeed, blackbody fits to the 
thermal spectra of old millisecond pul sars find inferred po- 
l ar cap radii of order 0.1-0.2 km dBogdanov et al.l [20061 ; 
iForestell et al.l 120141) . Using hydrogen atmosphere models, 
emitting rad ii of ~ 1 km, more typical of NS polar caps, are 
found (e .g. IZavli n et al.ll2002t iBecker & Aschenbachl120021 : 
iBogdanovet all 120061) . We therefore try fitting the thermal 
component with the NSATMOS model again, now permitting 
radiation from only a small fraction of the stellar surface, by 
freeing the normalization parameter. This result of this fit 
gave an aphysical normalization (corresponding to the frac¬ 
tion of the neutron star emitting) 2> 1. Thus, so fitting was re¬ 
peated with an alternate neutron star atmosphere NSA model 










































4 


Arnason et al. 


(IZavlin et al.ll 1996b . with magnetic field set to 0 G, which is a 
reasonable approximation for the 10 8-9 G regime (NSA + PL 
in Table [2]). The resulting fit is again an improvement over the 
single power-law fit (x 2 /dof = 187.99/188, F-statistic of 3.55 
and probability 3.1 x 10“ 2 of this improvement compared to 
the variable Nh single power-law model being due to chance). 
The best-fit temperature is 0.61 q 2 keV and the inferred radius 

of the emitting region is 0.3 /q/ km, with the temperature up¬ 
per error limit pegged to the upper boundary of the model. 

Finally, we tried a broken power-law model, where the 
spectrum is described by a power law of index F | = 1.37// 
up to an energy of 2.7+gg keV, where it becomes a power law 
of index T 2 = 1-9^0 2 - This model also gives a significantly 
better fit than the single power-law (A\ 2 = 9.89 over the 
single power-law, F-statistic 5.02, probability 7.5xl0~ 2 of 
this improvement compared to the variable N H single power- 
law model being due to chance). Releasing the power-law 
indices or break energies between observations did not im¬ 
prove this fit. The broken power-law model would presum¬ 
ably imply a bremsstrahlung origin of the emission, with 
the break energ y sugges ting the bremsstrahlung temperature 
(IChakrabartv et al.ll2014l) . so it is reasonable to search for a 
separate emission contribution from the neutron star surface 
in this case. However, as above, adding a NSATMOS compo¬ 
nent (with emitting size fixed to the whole neutron star) gave 
no improvement, and adding a BBODYRAD component gave 
very little improvement (A^ 2 = 0.01). Although it is an em¬ 
pirical fit only, we adopt the broken power-law model as our 
best-fit as it gave the best improvement compared to the vari¬ 
able Nh single power-law model. 

3. OPTICAL DATA ANALYSIS 
3.1. WFC3 Photometery 

The 2012 HST observation, taken in one orbit with the 
WFC3 camera using the UVIS channel in subarray mode, 
one day before a Chandra observation, consists of exposures 
in WFC3 B (F438W, total exposure length 1360s), broad V 
(F606W, total exposure length 188s) and I (F814W, total ex¬ 
posure length 332s). The exposures were obtained using the 
standard 4-point box dither pattern, with integer plus subpixel 
dithers, to allow forbad pixels and increase the effective angu¬ 
lar resolution. A lkx lk subarray was chosen to minimize the 
overhead time and thus allow the necessary exposures to fit 
into a single orbit. In programming the observations, care was 
taken to constrain the roll angle in order to prevent diffraction 
spikes from a nearby bright giant from falling on the M15 X- 
3 counterpart. The frames were combined with correction for 
cosmic rays using the IRAF/STSDAS Multidrizzle task with 
2 x oversampling, giving an effective pixel size of 0.02". Pho¬ 
tometry on the drizzle-combined frames was carried out using 
the DAOPHOT/ALLSTAR software package, since the field 
of M15 rapidly becomes crowded near the core. PSFs were 
constructed using 100 candidate PSF stars per filter. Photo¬ 
metric calibration to the VEGAMAG system was performed 
by doing aperture photometry on the PSF stars within a 0.1" 
diameter aperture, finding th e ape rture correction to an infinite 
radius aperture from lSirianni et all ( 2005 ). calculating the off¬ 
set between the ALLSTAR photometry and the aperture pho¬ 
tometry, and applying the calibrations from the HST WFC3 
calibration webpageQ 

6 http ://www. stsci. edu/hst/wfc3/phot_zp_lbn 
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Figure 2. X-ray spectra of M15 taken from the four ACIS-S Chandra obser¬ 
vations, fitted to a single power-law with column density as a free parameter. 
The spectra are sorted with the most recent at the top. Note that many of the 
residuals lie above the model between 1 - 3 keV and all of the residuals lie 
below the model above 3 keV, indicating the potential presence of curvature 
in the spectrum (particularly for the MJD = 56188 epoch). These data have 
been rebinned for clarity. 



Figure 3. Finding charts for the optical counterpart of Ml5 X-3. N is up and 
E is to the left. The field of each frame is 1.4" wide by 1.6" tall. In each 
frame, the location of M15 X-3 is identified with a red arrow. All images 
are drizzle-combined frames from the HST/WFC3 observation taken in 2012 
(MJD=56187), nearly simultaneous with Chandra ObsID 13710. 
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Table 2 

Summary of fits to Chandra observations of M15 X-3 


Model 

SPL 1 

SPL 2 

BB + PL 

DISKBB + PL 

NSATMOS + PL 

NSA + PL 

BPLf 

Nh 

[10 20 cm- 2 ] 

[4.6] 

9+ 

6+4 

+5 

10+ 

5+ 5 

J -5 

6+ 

r. 


1 42"*" 0 - 03 
^ -0.03 

i ^7+0.07 

' —0.07 

1-5+2 

1+0.5 

1 C7+0.04 

AJ ' -0.07 

1 5 +0 - 3 

1 ‘^—0.3 

1 3 +01 

1 0.2 

r 2 








l 9+ 02 

1 '^—0.2 

Fbreak 

[keV] 







9 7+0.4 
z,/ -0.6 

Chennai 

[keV] 



0 7+ 0 - 3 
u>/ -0.2 

A ' Z -0.3 

0 05” 1 ” 0 ' 03 
u ' u -0.04* 

0 6+ 0 - 3 * 
uo -0.2 


^thermal 

[km] 



0 2 -1-0 ' 1 
u - —0.1 

0 7”*” 0 ' 3 
y - / -0.5 

[10] 

0 3”*” 0 ' 4 
U.J-0 2 


L x , MJD = 55069 

[10 34 ergs” 1 ] 

n q'2+0.04 

-0.04 

O oci+0 04 

U '° - 0.04 

0 77+ 0 05 

u,// -0.05 

0 74+ 0 - 06 
u -'^-0.06 

O oa+0 04 

U.5J_oo4 

0-761“- 

0 75+ 0 05 
u - /:) -0.05 

PL fraction, MJD = 55069 


[1] 

[i] 

U -oo_0.09 

0 7 +0 - 2 
u/ -0.3 

0-99+;°o! 

0.841+ 4 

[i] 

L x , MJD = 55071 

[10 34 ergs” 1 ] 

0 00+0.0/ 
u - vy -0.07 

0-98++ 

n qt}+0.08 
0.08 

0-91+09 

0.98+;° 0 7 

0 92+ 0 - 08 
U ‘ yz -0.08 

n 90+ 0 08 
u -^ u _o.07 

PL fraction, MJD = 55071 


[1] 

[1] 

0 9o+ 0 06 
uyu -0.08 

0-8+2 

0 99+ 0 01 
u yy -0.01 

O 07+0-03 
U -o/_0.04 

in 

Lx, MJD = 55097 

[10 34 ergs -1 ] 

Q Q2 + 004 

^ -0.04 

Q Q2 + 0 04 

u -^ z — 0.04 

O Q7+0-05 
U ° 7 - 0.05 

O O In +0-06 

-0.07 

0 92+ 0 04 
u ‘^-0.04 

0-861“- 

0-84+;“ 

PL fraction, MJD = 55097 


[1] 

[1] 

0 89+ 0 07 
0.09 

0 7 +0 - 2 
u -0.2 

0 99+ 0 01 
u -^-o.oi 

0-861+j 

[i] 

Lx, MJD = 56188 

[10 34 ergs” 1 ] 

1 o +01 

A - u -o.i 

A - u -o.i 

A * u -0.1 

i-o+;! 

1 1 +0J 

Ai -0.1 

1 O” 1 ” 01 
A - u -o.i 

1 o +01 
i,u -0.1 

PL fraction, MJD = 56188 


[1] 

[1] 

0 91 +0.06 
u ■^ i _0.08 

0.8+-J 

0 99+ 0 01 
u - y ^-o.oi 

n rr+^'^ 3 

u -oo_0.04 

[i] 

X 2 /d.o.f 


214.13/191 

195.08/190 

187.52/188 

188.04/188 

195.03/189 

187.99/188 

185.19/188 

n.h.p. 


0.12 

0.38 

0.50 

0.49 

0.37 

0.49 

0.54 


Note. — M15 X-3 best spectral fit parameters, all ACIS-S data. Values in square brackets are fixed as part of the model choice. Error values marked with * 
indicate a model limit. In this table SPL 1 is the single power-law with column density fixed, while SPL 2 is the same but with column density free, f We adopt 
the broken power-law model as our best-fit. 



B-V 


V-I 


Figure 4. Color-magnitude diagram for Ml5, from 2012 //S77WFC3 data. 
The red dot indicates the location of the optical counterpart for M15 X-3, 
with errors calculated based on the half-width of the main sequence at the 
magnitude of the counterpart. Note that M15 X-3 appears to be bluer than 
the main sequence, especially as one moves to bluer filters. Aside from M15 
X-3’s optical counterpart, only stars with small formal photometric error (< 
0.025 mag) were kept to clean the image. 


As shown in Figure [3] the optical counterpart of M15 X-3 
is clearly detected in all 3 filters. Its magnitudes in this epoch, 
as well as in previous usable HST epochs, are given in Table 
[3 CMDs of the cluster were constructed in B - V and V-I, 
as shown in Figure 0] where the counterpart is bluer than the 
main sequence and becomes progressively bluer as one moves 
to bluer filters. 


3.2. PysynphotModelling 


Table 3 

Summary of HST Observations of M15 X-3 


Proposal ID 

Date 

MJD 

Filter 

Instrument 

Magnitude 

5742 

1994-10-26 

49651 

F336W 

WFPC2 

21.5T0.2 

9039 

2002-04-05 

52369 

F555W 

WFPC2 

22.0 ± 0.2 

9039 

2002-04-05 

52369 

F439W 

WFPC2 

23.7 ±0.8 

12751 

2012-09-17 

56187 

F814W 

WFC3 

21.69 ±0.04 

12751 

2012-09-17 

56187 

F606W 

WFC3 

22.34 ±0.09 

12751 

2012-09-17 

56187 

F438W 

WFC3 

22.77 ±0.12 


Note. — Observations in 2012 were calibrated to the VEGAMAG system, 
whi le observat ions in 1994 and 2002 were calibrated to the STMAG system, as 
per lHeinke et all 120091) 


To determine the nature of the M15 X-3 companion using 
its o ptica l _colors, the system is modelled using Pysynphot 
(iSTSel development Team! 12013 1. Previous observations of 
M15 X-3 in the X-ray and optical bands have suggested that 
the system is likely to be a neutron star with a low-mass main 
sequence companion dHeinke et al.ll2009l) . M1 5 X-3 appears 
to be too X-ray luminous to be a WD accretor (IVerbunt et alJ 
119971) and its X-ray/optic al flux ratio rules this out conclu¬ 
sively dHeinke et al120091) . The optical spectrum of the M15 
X-3 system is therefore modelled with two components: a 
low-mass main sequence companion and an accretion disk 
whose emission is approximated by a power law. 

The range of power-law indices used to model a power- 
law accretion disk (and its expected optical spectrum) de¬ 
pends on two limiting behaviours. For small frequencies, 
F v approaches a Rayleigh-Ieans description where F v oc v 2 . 
For a conventional disk, the spectrum of the outer part of 
the disk behaves approximately as F v oc v V 3 dFrank et all 
120021) . In the wavelength space utilized by Pysynphot for 
spectral modelling, this corresponds to a photon index Ta = 
—4 in the Rayleigh-Jeans limit, and Ta = —7 /3 for a con¬ 
ventional disk. In fitting the power law, I’a indices be¬ 
tween —4 and —2 were iconsidered. To have an index 
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larger than —2 (which corresponds to the peak of the spec¬ 
trum) would imply that the measured WFC3 bandpasses 
lie close to the Wien limit of the spectrum. Observations 
of quiescent disks in other accreting NS and black hole 
systems (e.g. iFroning et al.H201 It iHvnes & Robinsonl 120121 
iWang et al J 120 1 31) indicate that the disks have typical inner 
temperatures between 5000 and 13000 K. Compared to Aql 
X-l in quiescence, M15 X-3 is mo re luminous in both the X- 
rays (by a fa ctor of 20: iNar avan et a l .11 1 9971) and UV (by a 
factor of ~4; IHvnes & Robinsonll2012l) . Assuming they have 
a similar temperature-radius relation, we conclude that M15 
X-3 is likely hotter than Aql X-l. Since Aql X-l’s disk tem¬ 
perature was 12400 ± 1400 K (IHvnes & Robinsonil2012l) . it is 
unlikely that the disk temperature of M15 X-3 is low enough 
to place the Wien tail into our WFC3 bandpasses. 

The power-law component is processed with Synphot’s 
built in power-law model, while the low-mass companion is 
modelled using stellar models from the Castelli and Kurucz 
library of stellar atmospheres included in the Synphot pack¬ 
age. The advantage of this library is that it permits the mod¬ 
elling of main-sequence stars at the low metallicity charac¬ 
teristic of M15. As inputs, the Castelli and Kurucz models 
require the following: effective temperature T e s, logarithm 
of surface gravity log(g), and metallicity [M/H], To obtain 
these input parameters for a mai n-seq uence star of a particular 
mass, results were taken from lBaraffe et al.l (119971) model cal¬ 
culations for low-mass stars at [M/H] = —2.0, which ^rea¬ 
sonable given that M15 has [Fe/H] = —2.34 (ICarretta et al. 
12009 ). The nature of the companion and disk is determined 
by fitting the B, V, and I magnitudes predicted by (Synphot) 
synthetic photometry for a particular choice of low-mass com¬ 
panion and accretion disk against the measured WFC3 mag¬ 
nitudes for the M15 X-3 system using a standard chi-square 
test. A grid-search for the best chi-square value is performed 
over the parameter space for each free component — mass in 
the case of the companion, and power-law index in the case 
of the accretion disk. Four cases were considered: A low- 
mass main-sequence companion only, a power law accretion 
disk only (with free index), a companion and power-law disk 
with a free index, and a companion and power-law disk with 
index fixed to —7/3. In the case of a low-mass companion 
+ free index power-law disk, testing only at the mass values 
given in lBaraffe et al.l (119971) gives a best-fit companion mass 
of 0.5±0.1 Mq, with a power-law index Ta = —3.28. To 
estimate a more precise mass, linear interpolation of model 
parameters was performed on the table of masses in the range 
0.200-0.700 M q . 

As shown in Table [4] neither a power-law disk description 
nor a companion alone is sufficient to explain the observed op¬ 
tical colors of M15 X-3. The best fit is obtained by having a 
combined spectrum, where the accretion disk of M15 X-3 and 
its low-mass companion both contribute to the optical flux. 
Allowing the mass of the companion as well as both power 
law parameters (Ta and the normalization) to vary gives an 
acceptable fit with an inferred companion mass of 0.5 IS/)®™ 
Mq. However, with this many free parameters, the degrees of 
freedom reduce to zero, overconstraining the model and lim¬ 
iting its utility. Such a fit is only useful for constraining the 
model parameters — in this case, the fit shows the stellar com¬ 
ponent is not consistent with a low stellar mass. If the index is 
fixed at the conventional accretion disk value of —7/3 an ac¬ 
ceptable fit is also obtained, with an inferred companion mass 
of 0.440/QQgQ Mq. The inferred spectrum of this combined 



Figure 5. Inferred optical spectrum for M15 X-3. The spectrum of the 0.440 
Mq companion is plotted in orange. The spectrum of the accretion disk, 
modelled as a power-law, is plotted in green. The combined spectrum is 
plotted in cyan, with the black points representing measured optical fluxes in 
each of the three HST filters. The bandpasses of each of the three HST filters 
are plotted in arbitrary units in purple for reference. 


system is plotted in Figure 0 This choice of index describes 
the accretion disk as the largest contributor to the system’s 
optical flux in B and V, while in I the flux contribution of the 
disk and companion become approximately equal. This nat¬ 
urally explains why M15 X-3’s optical counterpart appears 
much bluer than the main sequence, especially in bluer fil¬ 
ters, as it is the accretion disk which dominates emission at 
bluer wavelengths. This fit has x 2 = 0.44 for 1 degree of free¬ 
dom. With this small d.o.f., this low ;y 2 corresponds to a null 
hypothesis probability of 0.49 (close to the null hypothesis 
probability of 0.5 for a reduced \ 2 = 1 for large d.o.f.), which 
indicates that there is no need to consider additional system¬ 
atic errors when constraining our fit. Note that the fits prefer 
a power-law photon index of —4, consistent with a relatively 
high temperature for the accretion disk. 

3.3. X-ray Irradiation of the Companion 

Given the relative blueness of the optical counterpart of 
M15 X-3 compared to other members of the cluster, it is rea¬ 
sonable to ask whether its optical colors may be substantially 
altered by X-ray irradiation of the main sequence companion 
of M15-3 by the primary. To determine whether X-ray irradi¬ 
ation is important, we compare the amount of X-ray light re¬ 
ceived by the companion star (measured near-simultaneously) 
with the total flux from the star. We begin by assuming a non- 
irradiated companion, M 2 = 0.44 Mq and a IAMq primary 
that radiates ss 10 34 erg s _I of X-rays. With a mass ratio 
q = AT/Afj one can write the binary separation as (follow¬ 
ing iFrank et ahll2002l) : 

a = 3.5 x 10 10 (M])5(1 +< 7 ) 5 (Phr) 4 cm. (1) 

For a lower MS star filling its Roche Lobe, the period can 
be approximated by 

M 2 Ri O.llPhr, (2) 

which means that the orbital perio d of the system is approx¬ 
imately 4 hours (Frank et al.ll2002l) . Substituting this relation 
in for M 2 gives the following: 
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Table 4 

Summary of fits to HST observations of M15 X-3 


Model 

r A 

U 1994 

B 2002 

#disk/#total 

V2002 #2012 

V 2012 

ho 12 

Inferred M co mp 

X 2 /dof 

Companion only 

PL only 

Fixed PL + companion 
Free PL + companion 

-2.001“;“* 

1-|] 

-4.ool 2 ;“: 

[0] 

HI 

0.98l“;“; 

0.94l“;“ 3 

[0] 

HI 

0 92+ 0 04 

0.04 

-0.03 

[0] 

in 

0.82l“;“ 7 7 

0.62l“; 34 

[0] 

[1] 

U ‘°-0.06 

0 71+0-28 
u - 11 -0.07 

[0] 

[1] 

0.68l“ : ; 3 

o.35l“;- 

[01 

[1] 

0.531-11 

O-lOlaoi 

0.540 ±0.005 

0 440+ 0 035 
u.44-u_ 0 060 

0 515+0.005 

U.JIJ_ O2 70 

33.78/2 

8.87/1 

0.44/1 

0.18/0 


Note. — Summary of fits to 2012 HST data of M15 X-3 optical counterpart using PYSYNPHOT modelling. Values in square brackets were fixed 
because of the model choice, while uncertainties marked with a single asterisk indicate that the uncertainty is at the hard limit of the model. Note 
that adding a power-law with free index to the model adds two model parameters (PL normalization, PL index), while a power law with a fixed index 
adds only one model parameter (PL normalization). A low-mass main-sequence companion adds only one free parameter to the model (mass of the 
companion). 


a ~ 3.5 x 10 lo Mi(l +q)* cm ' ^ 

For M15 X-3, this gives a~1.3x 10 11 cm. This means that g 
the X-ray flux received by the companion at its surface (as¬ 
suming isotropic emission) is 

F x « = 3.5 x 10 10 erg cm -2 s” 1 . (4) 

47ror 

By contrast, the total flux of X-3’s companion pan be esti¬ 
mated by taking a lumino sity of L+ r= 0.05L -., dBaraffe et al.l 
119971 ) and R+ « O.44^ 0 ([Frank etaTl 120021) to get the com¬ 
panion’s flux at its surface to be: 

F-k ~ T“ 77 r = ! -6 x 10 ' 0 erg cm -2 s -1 . (5) 

4nR + 

Since the received X-ray flux is roughly a factor of 2 larger 
than the flux emitted by the companion itself. X-ray irradia¬ 
tion is likely to be a contributing factor to the optical color 
of this system. This means that it is likely that the com¬ 
panion of M15 X-3 is less massive than the fitted mass of 
0.440^q[)6o Mq. To constrain the mass of the companion, we 
assumed that the companion is irradiated by its primary with 
Tprimaiy = 10 34 erg s^ 1 and re-radiates these X-rays away as a 
blackbody of temperature 

T ( 17 A 4 UL\ 

T irradiation — ( T2 ' 77~ I (6) 

\F*+F X J 

We assume that the emission is dominated by the irradiated 
half of the star, and thus that: 



Modelling the star as a blackbody with temperature T and 
radius R, the acceptability of various fits was explored over 
the parameter space in Figure[6] When compared with a max¬ 
imally irradiated main sequence star, at the 3cr level we can 
still rule out the possibility of a 0.15/Jq or smaller compan¬ 
ion. This rules out a brown dwarf, white dwarf, or planetary 
companion. At the 2a level, the radius of the companion lies 
between 0.19 and 0.41 R@, which suggests a ~ 2 — 4 hour 
orbital period. 



Radius (i? Q ) 

Figure 6. Contour plot of the fit of a single power law and irradiated com¬ 
panion (represented as a blackbody) for varying choices of the companion 
size and temperature. The orange, purple, and cyan curves indicate the 1 a, 
2<t, and 3<r curves, respectively. The bottom green curve indicates the tem- 
perature/radius relation for non-irradiated low-mass main-sequence stars of 
metallicity [ M/H ] = —2.0, while the top green curve represents the temper- 
ature/radius relation for the same stars maximally irradiated by 10 j4 erg s' 1 
of X-rays. For reference, the black point indicates the best-lit 0.44/V/q, star 
for the non-irradiated case. Regions in grey are excluded for physical reasons 
- those below the non-irradiated main sequence are aphysically colder than 
a main-sequence star of the same radius, while those above the maximally 
irradiated main sequence are hotter than possible for a main-sequence star 
irradiated by only a primary with an X-ray luminosity of 10 34 erg s — *. 


4. DISCUSSION 
4.1. Nature of System 

M15 X-3’s unusual X-ray behavior for a globular cluster 
X-ray source, and position blue-ward of the main sequence, 
suggests that we check whether it can be explained as a fore¬ 
ground or background source. If M15 X-3 were a foreground 
star, we would expect it to be redder, rather than bluer, than 
the main s equence. We use the empirical relations from 
iBaldi et akl (120021) to estimate the probability of finding a 
background AGN at or above the flux of M15 X-3 projected 
within the half-mass radius of M15. For M15 X-3, we use the 
bright state flux of S ~ 5.6 x 10“ 13 erg/(s cm 2 ), which gives an 
expected 0.6 sources per square degree, or 5 x 10~ 4 sources 
expected in the half-mass radius of M15, roughly the region 
where M15 X-3 was discovered. (Note that the likelihood of 
finding a source at the distance of M15 X-3 from the cluster 
center is a factor of 9 smaller.) This suggests that M15 X-3 is 
unlikely to be a foreground or background source. 
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Figure 7. Long term light curve for M15 X-3. ROSAT, Chandra HETGS, 
and Chandra HRC-I data (where M15 X-3 is in quiescence) are taken from 
IHeinke et alJ <20091) . while the Chandra HRC-I (where M15 X-3 is bright) 
and ACIS-S data were analyzed in this paper. The vertical dashed lines in¬ 
dicate HST observing epochs, including the 2012 observation that is nearly 
simultaneous with an ACIS-S observation. 


We can quantitatively consider the probability of the opti¬ 
cal counterpart being an unaffiliated interloper by consider¬ 
ing the local density of optical sources that could be potential 
counterparts. We consider a region of 5” around M15 X-3, 
as there is a strong density gradient of sources with distance 
from the center of the cluster. Within this region, there are 
3 objects that have V magnitudes in the range 19 < V <23 
that lie as far to the blue side of M15’s main sequence as M15 
X-3’s proposed counterpart, including the counterpart itself. 
This gives a chance probability of one of these objects falling 
within the M15 X-3 error circle (0.05" — the offset between 
the Chandra and HST relative astrometry, determined using 
the optical and X-ray detections of the nearby XRB AC 211 
IHeinke et al.ll2009i) of 3 * (0.05”/5”) 2 ~3x 10” 4 . Standard 
counting errors would give 3^ 2 g objects at la, and 3^'8 at 
3a. At the upper end, the probability of a non-associated op¬ 
tical counterpart that is a member of the cluster is f=s 0.13% 
even at 3a. The remaining possibility for the optical emission 
of M15 X-3 is that the reddish companion is a cluster mem¬ 
ber unassociated with M15 X-3 itself. However, if the B flux 
is dominated by M15 X-3, and the I flux is dominated by an 
unassociated main sequence star, the offset between these two 
objects could not be more than about 1 WFC3 pixel = 0.04", 
as the counterpart is well-aligned in the B and I frames. This 
gives a chance probability of an unassociated red star happen¬ 
ing to align with M15 X-3’s blue counterpart of « 3%, consid¬ 
ering stars within one magnitude of the proposed counterpart 
in I. 

As discussed in IHeinke et al.l ( 120091) , we can immediately 
rule out a white dwarf accretor on the basis of long-term 
variability and the X-ray to optical flux ratio. Radio ob¬ 
servations of M15 X-3 argue against the idea that we view 
it edge-on, permitting most of its intrinsic luminosity to be 
obscured by the accretion disk, because the radio emission, 
which scales wi th the intrinsic X-r a y luminosity, should still 
be visible (e.g.. iBower et al.l 120051) . IStrader et al.l (12012l) re¬ 
port no source at M15 X-3’s location, with an rms sensitivity 
of 2.1 /iJy/beam. A non-detection implies that M15 X-3’s ra¬ 


dio brightness cannot be higher than 10.5 /jJy at 5a. For a 
black hole at 10 36 erg s -1 , with X-rays obscured by a high 
inclination accretion disk, we would expect a radio flux of 
roughly 10 30 erg s_'. At 10.3 kpc , this would give roughly 
1 mJy at 6 GHz dGallo et al .1120121) . For a neutron star, we 
would expect closer to 20 p Jy at 6 GHz dMigliari & Fender! 

l200fih . 


4.2. X-ray Spectra 

M15 X-3 appears to be the faintest VFXT for which we 
have an X-ray spec trum of very high quality; fo r instance, the 
targets discussed in I Armas Padilla et al.l (I2013al) have Lx (0.5- 
10 keV) of 5-10xl0 34 erg s -1 , vs. ~ 1 x 10 34 erg s -1 here. 
The X-ray spectrum is somewhat atypical, and can be well-fit 
by either a broken power-law or a power-law plus a thermal 
component. The thermal component, if fit by a blackbody, 
produces an inferred radius of 0.2 +q'J km (or 0.3 ^q 4 km if fit 
with an NSA model), which is significantly smaller than the 
inferred blackbody radii (~5 km) of the thermal compo nents 
in the persistent VFXTs discussed by lArmaTP adilla et al. 
(l2013al) . It is unclear if the possible thermal component here 
has the same physical origin as in the other transients. 

The X-ray spectrum of material falling onto a neutron star 
is predicted to include a thermal component (from the lower 
atmosphere) plus a very hot surface layer that effectively 
Comptonizes the outgoing radiation (iDeufel et al.l200ll) . This 
model may provide a first step toward understanding the spec¬ 
tra of sev eral observed neutron star transients at low luminos¬ 
ity (e.g., lArmas Padilla et aI1l2013aL iBahramian et al.lf20l4 
iChakrabartv et alJ 120141) . In the case of M15 X-3, the to¬ 
tal luminosity and small inferred radius of the thermal com¬ 
ponent are consistent with such a high proton energy flux 
(~ 10 24 ergcm~ 2 s _1 ) that the Deufel et all (120011) model pre¬ 
dicts an entirely non-thermal spectrum. To reproduce the ap¬ 
proximately 10% thermal component in M15 X-3, one might 
invoke a central high proton flux hotspot that is smaller than 
the radius of the thermal component we measure, where the 
thermal component arises from the surrounding area being 
heated either by conduction or by a lower proton energy flux. 
It would be of great interest to obtain models such as Deufel’s 
in a convenient form for the spectral fitting of X-ray transients 
0 . 


4.3. Explanation for VFXT Behavior 

If M15 X-3 spends most of its time in its "bright" state, with 
a persistent luminosity of 10 34 erg s _1 , then we can infer an 
upper limit on the time-averaged accretion rate by assuming 
that all of its 0.5-10 keV X-ray luminosity is from accretion. 
This assumption implies that M15 X-3’s time-averaged accre¬ 
tion rate could be at most 9 x 10~ 13 Mq/ year. Under conven¬ 
tional binary evolution, the low X-ray luminosity of M15 X-3 
is difficult to explain. One suggestion is that the low accre¬ 
tion rates needed for VFXTs is possible from a brown dwarf 
or planetary-mass companion. Another option is a system 
where the accretor is an intermediate-mass (~1000 Mq) black 
hole with a pri mordial companion that i s very low mass in the 
current epoch (King & Wiinandsll2006l) . If the companion of 
M15 X-3 is a brown dwarf or other very low-mass companion, 
then we would expect the optical emission from the system to 
be produced only by the accretion disk. At 10.3 kpc, a brown 

7 The zamp models o UZampieri et all 119951) do not include the hot Comp- 
tonizing layer of the atmosphere. 
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dwarf companion would be well below the detection thresh¬ 
old from the HST observations, in contrast with the presence 
of a second, non-accretion disk component. The poor fit of a 
power law alone to the BVI values for M15 X-3 suggests that 
a disk-only model is insufficient to explain the nature of M15 
X-3 — the accretor’s companion is not sufficiently low-mass 
for this explanation to be tenable. The location of M15 X-3, 
21 " from the cluster center, places it well outside the region 
where a n intermediate-mass black hole would be e xpected to 
wander dChatteriee et al.ll2002tlGerssen et al]l2002l) . 

It is also plausible for some VFXTs that accretion from the 
wind of a MS star could explain low accretion rates. In the 
case of M15 X-3, a companion of M < 0.44M is unlikely 
to be driving wind-fed accretion i n this system (IHeinke et all 
120141) . It has also been argued bv lin’t Zand et all (120071) that 
persistently low accretion rates require a very small, com¬ 
pletely ionized disk, implying an ultracompact system. Our 
optical observations of M15 X-3 are inconsistent with an ul¬ 
tracompact system, as we have found evidence for a relatively 
large (0.15-0.45 Rq) donor. 

An alternative explanation for the unusually low X-ray per¬ 
sistent X-ray luminosity of M15 X-3 could be that it is an 
otherwise normal LMXB (with a 2-4 hour orbital period) 
where the NS magnetic field is sufficiently strong to inter¬ 
fere with accretion, via the propeller effect (in which the rota¬ 
tion of_dieNS^sjnugneticfieldis able to repel inflowing mat- 
ter; llllarionov & Sunvaevll 19751) . ID’Angelo & Snruiti d2012l) 
modelled a leaky propeller and "trapped" disk, leading to con¬ 
tinuous low-level a ccretion; similarly leak y propellers have 
been s t udied by e.g. iRomanoya eLal. (2003),^st^ugovaetalJ 
(120061) . iKulkarni & Romanoval (12008 ). IHeinke et all (120141) 
argued that all the quasi-persistent VFXTs such as M15 X-3 
are in this physical state. This argument is based, on one hand, 
on the difficulty of explaining continued extremely low-level 
accretion via disk instability models, and on the other hand 
the example of transitional millisecond pulsars, which main¬ 
tain low-level accretion and clearly show evidence of mag¬ 
netic effects on the accretion flow. 

Three transitional millisecond pulsars (which cycle be¬ 
tween states of low-level accretion at 10 33 — 10 34 erg s 1 , 
non-accreting radio pulsar states at ~ 10 31 erg s _1 , and oc 
casio nal higher-Ly accreting states) have recently been iden 


casional nigiier -Lx accreting states) nave recently been iden- 
tified (|Ard]ibald_eialJ2QQ9t^|PaniUoeta]J2()13t iPatnmoetal 


2014: fStappers et al.ll2014UBassa et al.ll2014UBogdanov et al 


2014 ; iRov et al .1 120 1 4l) . The recent detection of X-ray pul¬ 
sations from two of these transitional millisecond pulsars in 
their Lx ~ 10 33 — 10 34 erg s” 1 state (lArchibald et all 12014! 
iPapitto et al.l2014l) strongly indic ates that these pulsars a re in¬ 
deed in this leaky propeller state (ID’Angelo et al.ll2014l) . The 
quasi-persistent VFXTs, including M15 X-3 , may therefore 
alsobe transitional millisecond pulsars dDegenaar et al.ll2014l 
IHeinke et al.111201 4) . At least one transitional millisecond pul¬ 
sar has been o bserv ed tobe relatively radio-bright in its X-ray 
active phase (iDeller et al .1120141) . Of the known tMSP sys¬ 
tems, M15 X-3’s X-ray luminosity is closest to XSS J12270- 
4859, which would imply it has a radio brightness of roughly 
5 x 10 27 erg s -1 at 5 GHz if it is a similar source. Assuming 
a spectrally flat source, M15 X-3 would have a 5 GHz flux 
density of approximately 5.3 pJy at 10.3 kpc, roughly a factor 
of two too faint to be detected by existing radio observations. 

A relatively conservative estimate of M15 X-3’s duty cy¬ 
cle is roughly 70%, as seen in Figure [7] However, the quasi- 
persistent VFXTs with a duty cycle resembling that of M15 


X-3 tend to have relatively soft power-law fits, or harder 
power-law fits only when a t hermal component is also present 
(1 Armas Padilla et af1l20 1 3al) . The X-ray spectrum of M15 X- 
3 is unlikely to be an intrinsically harder spectrum being sup¬ 
pressed at lower energies because of absorption effects — ra¬ 
dio observations suggest the system is not being viewed edge- 
on, and the column density of M15 X-3 we measure is only 
< 10 21 cm -2 . 


5. CONCLUSION 

We have performed a combined optical/X-ray followup ob¬ 
servation of the VFXT M15 X-3. Our X-ray observation, 
and other archival Chandra observations, found M15 X-3 in 
its "bright" state, persistently at a luminosity of « 1 x 10 34 
erg s -1 . The X-ray spectrum, if fit with a power-law, show 
strong residuals. The X-ray spectra are better fit by a power 
law of index 1.28 +q^ up to an apparent spectral break at 
2.7ly g keV, where it is described by a power law of index 
1 -9^0 2 - A power-law plus thermal component also describes 
the spectrum acceptably, giving an emitting region of ~ 0.21 
km for kT ~ 0.6 keV (for a blackbody model) or ~ 0.4 km 
for kT ~ 0.6 keV (for a NS atmosphere model). The HST ob¬ 
servation clearly detects a blue optical counterpart, which in¬ 
cludes contributions from both the system’s accretion disk and 
a companion with emission consistent with a main-sequence 
star of mass 0.440^'^ Mq. The actual mass may be lower, 
in the range 0.20-0.45 M 0 (at 2cr), since X-ray irradiation of 
the companion can cause the star to emit more strongly in 
the optical bands than it would as an isolated main-sequence 
star. M15 X-3’s persistently low accretion rate is difficult to 
explain in conventional binary evolution models. Combining 
this low accretion rate, with the evidence of the donor na¬ 
ture presented here, argues for accretion via a leaky propeller 
mechanism, as has been recently suggested for the transitional 
radio pulsars. 
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